Electric dipole-free interaction of visible light with silver metadimers 
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In subwavelength-sized particles, light-induced multipole moments of orders higher than the elec- 
tric dipole are usually negligibly small, which allows for the light-matter interaction to be accurately 
treated within the electric dipole approximation. In this work we show that in a specially designed 
meta-atom, a disc metadimer, the electric quadrupole and magnetic dipole can be the only excitable 
multipoles. This condition is achieved in a narrow but tunable spectral range of visible light both 
for individual metadimers and for a periodic array of such particles. The electromagnetic fields scat- 
tered by the metadimers fundamentally differ from those created by electric dipoles. A metamaterial 
composed of such metadimers will therefore exhibit unusual optical properties. 



In natural materials, effects of light-matter interac- 
tion, such as optical refraction and reflection, are nearly 
completely governed by the electric dipole excitations in 
the atoms and molecules composing the medium. While 
higher-order multipoles, such as electric quadrupoles and 
magnetic dipoles, can give rise to various optical phenom- 
ena their contribution to the overall interaction is 
very small. In optical metamaterials, however, light can 
excite also significant magnetic dipole moments, which 
can lead to such extraordinary optical phenomena as neg- 
ative refraction and near-field focusing 

Recently, it was shown that one can efhciently excite 
both magnetic dipole and electric quadrupole moments 
at visible frequencies in a pair of silver bars Q . In spite 
of the fact that the bars had rather large dimensions, 
the electric dipole contribution to the light scattering 
was still dominating in this structure. A metamaterial, 
in which the meta-atomic electric dipole moments would 
have a negligible contribution to the light scattering com- 
pared to the higher-order multipoles would be a unique 
optical material that has been neither observed in nature 
nor studied in scientific experiments. The radiation pat- 
terns of electric quadrupoles and magnetic dipoles funda- 
mentally differ from those of electric dipoles. Therefore, 
light propagation phenomena in such higher-order meta- 
materials must be treated in a conceptually different way 
than in an ordinary medium. As an example, according 
to the recently introduced electromagnetic boundary con- 
ditions in the presence of electric quadrupoles, the tan- 
gential components of the electric and magnetic fields can 
be discontinuous across the material boundaries [3, 
Furthermore, if the electric polarization P in the mate- 
rial could be made negligible, the interaction of light with 
the material would be described by Maxwell's equations 
in which P is absent. Such a material can exhibit reduced 
optical reflection and refraction at the interface with vac- 
uum. If, in addition, the absorption could be made low, 
the material itself could become essentially invisible. It 
is also remarkable that, in contrast to electric dipoles, 
the fields radiated by light-induced electric quadrupoles 
and magnetic dipoles in the forward and backward di- 
rections oscillate out of phase. This property could be 



used for creation of ultrathin optical elements, such as 
beam splitters and phase retarders, with unusual optical 
characteristics. 

In this Letter, we propose a concept of metamaterial 
that interacts with light solely via electric quadrupole 
and magnetic dipole excitations. These electromagnetic 
multipoles are of the same order, and for brevity we refer 
to this type of material as quadrupole metamaterial. To 
facilitate a practical realization of this concept, we search 
for a unit structure for the material, a meta-atom, that is 
(i) simple, such that the material can be relatively easily 
fabricated, (ii) much smaller than the visible-light wave- 
length, for the metamaterial to be treatable as homoge- 
neous in the visible spectral range, and (iii) symmetric, 
to enable a polarization independent response of the ma- 
terial to the field at normal incidence. To reach this 
goal, we employ a metadimer consisting of two closely 
spaced metal nanoparticles. Such dimers have previously 
been studied for excitation of antisymmetric [6|-|9| or dark 
plasmon modes that are characterized by simul- 



taneous oscillation of two oppositely directed electric cur- 
rents. The currents in those dimers, however, were only 
partially opposite and did not have equal effective am- 
plitudes. As a result, the dominating scattering by the 
structures still originated from electric dipoles. In order 
to exactly identify and evaluate the light-induced multi- 
pole moments in our meta-atoms, we use the electromag- 
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netic multipole expansion introduced, e.g., in Ref. 

The metadimer we propose as a unit of a quadrupole 
metamaterial consists of two axis-aligned silver discs, 
arranged as depicted in Fig. [T] The discs have the 
same thickness hi — h2 ^ 10 nm, but different radii, 
Ri = 15 nm and i?2 = 20 nm. The surface-to-surface 
separation of the discs is set to s = 10 nm. In our cal- 
culations, we choose the coordinate system to have the 
origin located on the dimer axis, at the center of the gap 
between the discs. The z-axis is directed such that the 
smaller disc is in the z < half-space, while the larger 
one is in the z > half-space. In order to have a scenario 
that can be realized in practice, the discs are considered 
to be embedded in a homogeneous dielectric of refractive 
index 1.5. This choice implies no loss of generality and 
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the conclusions to be drawn are equally applicable to any 
other choice of lossless surrounding. 
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FIG. 1. Illustration of the metadimer geometry. 

We assume that the metadimer is illuminated by a lin- 
early x-polarized plane wave propagating in the +z di- 
rection. Because of the chosen size and geometry for 
the metadimer, all moments of higher order than the 
electric quadrupole and magnetic dipole moments can 
be neglected. The symmetry of the metadimer with re- 
spect to the illumination direction and polarization en- 
sures that the electric dipole, magnetic dipole and elec- 
tric quadrupole moments are p = 'X-Px, = ym^, and 
= (xz -I- zx)g2;z, respectively. Here xz and zx are the 
outer products of the unit vectors x and z. 

We first calculate the scattered electromagnetic field 
distribution around the metadimer using the computer 
software COMSOL Multiphysics. The values of the elec- 
tric permittivity of silver are taken from Ref. We 
then expand the scattered field by using the multipole ex- 
pansion [13, [H]. The expansion coefficients are used to 
calculate the contributions C| , and Cj of the electric 
dipole, magnetic dipole, and electric quadrupole, respec- 
tively, to the scattering cross section of the metadimer. 
These modal cross sections describe the amount of opti- 
cal power radiated to the far-field by the corresponding 
multipole moment, relative to the intensity of the inci- 
dent plane wave. Thus, the modal cross sections enable 
us to compare the different multipole excitations. 

The modal cross sections of the metadimer as a func- 
tion of the vacuum wavelength Aq are depicted in Fig. [21 
Figure [2^ shows the electric dipole contribution to the 
scattering cross section over the whole visible spectrum. 
It can be seen that at around Aq = 594 nm, the electric 
dipole moment is suppressed and, as Fig. [2|d indicates, 
the power scattered by the electric dipole mode is neg- 
ligibly small compared to that scattered by the electric 
quadrupole and magnetic dipole modes. While the two 
remaining cross-sections C™ and Cj are small, they en- 
tirely determine the light-metadimer interaction at this 
wavelength. 

To obtain an intuitive picture of the electric current 
distribution in the metadimer, we calculate the electric 
dipole moments excited in each of the coupled discs. The 
magnitudes and phases of these moments are depicted in 
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FIG. 2. Spectra of the modal scattering cross sections of a sil- 
ver metadimer. The metadimer is embedded in glass and has 
the dimensions {R\,h\,R2,h2,s) = (15,10,20,10,10), in nm. (a) 
The whole visible spectrum of the cross section CJ due to the 
electric dipole mode, (b) The cross sections originating from 
each of the three lowest order modes around the wavelength 
of the electric dipole suppression. 



Fig. [3] In the spectral region around 594 nm, the electric 
dipole moments of the two discs oscillate out of phase 
with respect to each other (see Fig. [3^). A complete 
electric dipole suppression is obtained if, in addition, the 
two dipole moments have equal amplitudes. Figure [5)3 
shows that the magnitudes of the dipole moments are 
indeed equal at the wavelength of 594 nm, which is in 
agreement with Fig. [21 

We have verified by calculations that increasing the 
separation between the two discs will blue-shift the spec- 
tral location of the electric dipole suppression. For exam- 
ple, at a separation of s = 30 nm this location is shifted 
to a vacuum wavelength of 546 nm. Thus, varying the 
separation allows us to tune the spectral location of the 
dominating electric quadrupole and magnetic dipole scat- 
tering. 

Next, we show that in the multipole expansion of the 
electric current in the metadimer, the electric quadrupole 
and the magnetic dipole are interconnected excitations of 
the same order. The small size of the metadimer allows 
us to treat it as a point particle. Neglecting all multipole 
moments of higher order than the electric quadrupole and 
magnetic dipole, the electric current density of a point 
particle can be written in the form (see, e.g., [16]) 



J(r) = —iuj(p — ^ 



V)(5(r) -m X V(5(r). (1) 



At Aq = 594 nm, the dominating excitation in the 
metadimer can be seen as a second order electric current 
excitation, a current quadrupole. The major contribu- 
tion to the excited current quadrupole can be considered 
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FIG. 3. Spectra of (a) the phase and (b) the absolute value 
of the i-components of the excited dipole moments in the 
individual discs of the metadimer, normalized to the electric 
field amplitude Eq of the incident light. 



to assume the form 17 1 
EUr) 



U! 47re r 

(4) 

In contrast to an electric dipole, which radiates symmet- 
rically around its axis, the current quadrupole radiates 
primarily in the +2, and — z directions. Furthermore, the 
fields radiated in the +z and — z directions oscillate out 
of phase. The numerically calculated far-field radiation 
pattern of the metadimer, at the wavelength of the elec- 
tric dipole suppression, is depicted in Fig. For com- 
parison. Fig. shows the far-field radiation pattern of 
an electric dipole. The field profile in Fig. is in good 
agreement with Eq. ((4]). 





to consist of two opposite x-oriented current elements of 
length L separated by a distance s along the z-axis. If 
the element located at z > carries a complex-valued 
current of +Ixz and the one at z < a current of —Ixzi 
the electric current density in the point particle limit is 



Jxz{r) = -yiIxzLs—5[r). 

az 



(2) 



Comparing Eq. ^ with Eq. ©, we note that in the 
Jxz excitation the non-zero multipole moments are ^ = 
{5ci + i5t)iIxzLs I {2u}) and m = yIxzLs/2. Another con- 
tribution to the electric current excitation present in the 
metadimer is Jzxi in which we have a similar ^ = 
(xz -I- zx)z/2a;Xs/(2a;), but opposite m = —yIzxLs/2, 
where, however, Izx is small compared to Ixz- 

It is of interest to compare the field created by the cur- 
rent quadrupole excitation in the metadimer to the field 
created by an electric dipole. In spherical coordinates, 
the electric far-field of a point electric dipole oriented 
along the x-axis and situated at the origin of the coordi- 
nate system is 
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■ (0 cos cos (/) — sin 0) . (3) 



On the other hand, the electric far-field of the point cur- 
rent quadrupole introduced in Eq. ^ can be calculated 



FIG. 4. (a) Far-field radiation pattern of the metadimer at 
the wavelength of 594 nm. The excitation field is x-polarized. 
(b) Far-field radiation pattern of an x-oriented electric dipole. 
Arrows: electric field vectors. Colors: normalized intensity. 

To create a quadrupole metamaterial that can be 
treated as homogeneous, the metadimers should be ar- 
ranged in a lattice with a short period A << Aq. Thus, 
we have studied the metadimers arranged in an infinite 
two-dimensional periodic lattice. The electromagnetic 
coupling between the adjacent metadimers was found not 
to hinder the realization of the electric dipole suppres- 
sion, but simply to red-shift its spectral location. For 
example, for a square lattice with a period of 50 nm, 
we find the electric dipole suppression to take place at 
Aq — 618 nm. The electric current distribution (includ- 
ing both the conduction and polarization currents in the 
silver discs and in the surrounding dielectric) within a sin- 
gle unit cell of the lattice is depicted in Fig. [SJ The plot 
corresponds to the instant of time at which the electric 
currents have maximum values. The excitation in the 
metadimer exhibits zero electric dipole moment and is 
characterized by electric quadrupole and magnetic dipole 
moments of considerable strengths. In the figure, we also 
plot the amplitude of the magnetic field normalized to 
that of the incident wave. Within the metadimer, the 
magnetic field is seen to be significantly enhanced due to 
the strong opposite electric currents. These calculations 
clearly indicate the possibility to obtain a metamaterial 
with dominating electric quadrupole and magnetic dipole 



4 



excitations. 
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for the metamaterial layer to be treated as optically ho- 
mogeneous. 

We anticipate that practical realizations of the pro- 
posed metamaterial will not only help verify theoretical 
predictions of light interaction with higher-order multi- 
pole materials [3, [2^ , but also lead to discovering new 
optical phenomena. Furthermore, we believe that the ef- 
fects related to higher-order multipoles, such as optical 
activity, gyrotropic birefringence and dynamic magneto- 
electric effect , can be substantially enhanced and con- 
trolled in quadrupole metamaterials designed specifically 
for these purposes. 

This work was funded by the Academy of Finland 
(project 134029). 



FIG. 5. The electric current distribution and the magnetic 
field amplitude in an individual silver metadimer belonging 
to a two-dimensional square lattice of dimers embedded in 
glass with a period A = 50 nm. The illumination wavelength 
is Ao = 618 nm. Arrows: electric current density vectors. Col- 
ors: absolute value of the complex magnetic field normalized 
to that of the incident wave. 



As is evident from the symmetry of the metadimers, 
the response of the metamaterial layer is insensitive to 
the polarization of light at normal incidence. Another 
interesting property of the metadimer lattice is that the 
reversal of the illumination direction makes the dipole 
suppression phenomenon disappear and the lattice acts 
essentially as a lattice of electric dipoles. This can be 
used to create ultrathin bifacial optical components [3] ■ 

A three-dimensional quadrupole metamaterial can be 
obtained by stacking several layers of metadimers on top 
of each other. Owing to the geometrical simplicity of 
the proposed metadimers, the metamaterials composed 
of them should be relatively easy to fabricate. This would 
open up a possibility to study the properties of the still 
unexplored higher-order metamaterials also experimen- 
tally. 

In summary, we have introduced a concept of 
quadrupole metamaterials in which electric quadrupoles 
and magnetic dipoles entirely determine the light-matter 
interaction phenomena. We have proposed a realistic 
meta-atom, in the form of a metal-disc nanodimer that 
can be used to construct such a metamaterial for oper- 
ation at certain wavelengths of visible light. The spec- 
tral location at which the electric dipole contribution is 
suppressed can be changed by tuning, e.g., the disc sep- 
aration. Our calculations show that dominating higher- 
order multipole scattering by the metadimers can also be 
achieved in a periodic two-dimensional array. The size of 
the metadimers and the period of the lattice were chosen 
to be much smaller than the wavelength, which allows 
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